Critical Current of the Spin- Triplet Superconducting Phase in Sr2Ru04 
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There have been two different proposals for the spin-triplet order parameter of the superconducting 
phase in Sr2Ru04; an /-wave order parameter and the multigap model where some of the bands 
have the line node. In an effort to propose an experiment that can distinguish two cases, we study 
the behavior of the supercurrent and compute the critical current for these order parameters when 
the sample is a thin film with the thickness d ^ where ^ is the coherence length. It is found 
that the supercurrent behaves very differently in two models. This will serve as a sharp test for the 
identification of the correct order parameter. 
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Introduction: The order parameter of the supercon- 
ducting phase in Sr2Ru04 has been a subject of intensive 
research since its discovery in 1994 In particular, 

the field has been stimulated by the prospect of identify- 
ing the simplest electronic version of spin-triplet super- 
conductor. The early proposal made by Rice and Sigrist 
suggested the following spin-triplet order parameter with 
p-wave symmetry 3, 4J. 



A(k) = Ad {kx ± iky), 



(1) 



where A is the magnitude of the order parameter and d is 
a unit vector perpendicular to the spin of the condensed 
pair 1^ . This order parameter breaks time-reversal sym- 
metry and has a full gap on the Fermi surface. Indeed the 
flat ^^O Knight shift across 0, for the magnetic field 
along the ab plane, and the observation of spontaneous 
magnetic moment in /iSR are consistent with the spin- 
triplet order parameter and broken time-reversal symme- 
try, respectively. This encouraged theoretical investiga- 
tions of the effects of the characteristic collective modes 
(the spin waves and clapping mode) and the topological 
defects with the order parameter of Eq. 1 d, 0, 0; Ell • 

The experiments performed later on cleaner samples, 
however, reveal the existence of the nodal structure in the 
order parameter; the dependence of the specific heat 
[l^ , the T-linear behavior of the superfluid density flT 



peril 

the dependence of 1/Ti in NMR |ii|, the r2_behavior 



of the ultrasonic attenuation ll5|, and the V H depen- 
dence of the specific heat in a magnetic field ^3 ^^w 
temperatures. While these behaviors are consistent with 
an /-wave superconductor with nodes [l^ ItR E(Tj | . 

this may be at odd with the fact that the quasi-two 
dimensional system with strong ferromagnetic fluctua- 
tions would typically favor the p-wave superconductor 
with Eq.l I21L li^. In an effort to resolve this issue, the 
following order parameter with the horizontal node was 
proposed as a strong candidate JJj ■ 

A(k) = Ad (k^ ± iky) cos{ck^), (2) 

where c is the lattice constant in the c-axis. It was shown 



that the magneto-thermal conductivity would provide 
useful information about the validity of this order pa- 
ramete r l23[ |3| • Such experiments were indeed carried 
out |25ll26| . and the results are consistent with Eq.2, even 
though the data cannot exclude the possible presence of 
small amount of p-wave mixture. 

An alternative model was suggested by Zhitomirsky 
and Rice [13, . Their multigap model utilizes the fact 
that there exist three bands; a, /3, and 7 Q. Here the 
dominant p-wave order parameter with Eq.l resides on 
the active 7 band while the proximity effect leads to the 
following /-wave-like (it will be called /'-wave hereafter) 
order parameter with the line nodes in the a and P bands 



A(k) = Ad {kx ± iky) cofi{ckz/2). 



(3) 



The observation of a double transition in a recent specific 
heat measurement near Hc2 was interpreted in terms of 
this multigap model 29]. On the other hand, the specific 
heat and the magnetic penetration depth for low temper- 
atures (r ^ Tc) and low field {H ^ -ffc2) appear to be 
consistent with the /-wave order parameter of Eq.2. 

Thus it is not clear at the moment which order pa- 
rameter is the correct description of the superconducting 
phase; this is the fundamental issue for Sr2Ru04. Previ- 
ously it was suggested that the angle dependence of the 
magneto-thermal conductivity can be used to distinguish 
different order parameters 30] . The Raman spectra can 
be also used to detect different contributions from the 
clapping mode in the /-wave and the multigap models 
|3ll l32l |. These proposals, however, reply on the quanti- 
tative, albeit detectable, difference in two models. Thus 
it is still desirable to have a sharp test that leads to the 
qualitatively different outcome in two cases. 

In this paper, we provide such a proposal; the study of 
the critical current in two models for the spin-triplet su- 
perconductor. In particular, we investigated the behav- 
ior of the supercurrent for the /-wave and /'-wave order 
parameters. It is found that the supercurrent in the ab- 
plane has the same form in both cases and has significant 



2 



contributions from the quasiparticles. On the other hand, 
if the current flows along the c-axis, the supercurrent in 
the /'-wave is relatively immune to the quasiparticles at 
small current while that of the /-wave still acquires a 
large contribution from the quasiparticles. Notice that 
the supercurrent in the multigap model would be dom- 
inated by the p-wave component on the 7-band and its 
behavior would be the same as the case of s-wave super- 
conductor where the quasiparticle contribution is almost 
absent at small current [23 ■ Combining all these, it can 
be seen that the supercurrent in the /-wave case is much 
more affected by the quasiparticles. Thus the study of 
the supercurrent provides a clear mean to distinguish two 
different proposals for the order parameter of the super- 
conducting phase in Sr2Ru04. 

More specifically, we will consider a thin film with the 
sample thickness d ^ ^, where ^ is the coherent length. 
Under this condition, the magnitude of the supercon- 
ducting gap. A, and the supercurrent will be uniform 
across the system j^^l- When a uniform current flows, 
the Cooper-pair acquires a center of mass momentum 
of q^. If there were no contribution from the quasipar- 
ticles, the supercurrent would be simply proportional to 
Qs- The quasiparticles, however, change this behavior via 
the shifted quasiparticle dispersion, i?k, in the presence 
of the superflow 



(4) 



where = Vsk+^k ^^"^ "^k = 9^k/9k is the quasi- 
particle velocity. The total current, js, as a function of 
s = vpCls rises linearly at small s and has a maximum 
at some value of s, then drops as s is further increased. 
When djs/ds < 0, the system is unstable; thus the critical 
current is determined by the value of s where djs/ds = 
[3^ . The detailed behavior of the supercurrent and the 
value of the critical current in each case are discussed 
below. 

/-wave order parameter: We will consider two dif- 
ferent cases; the current in the a6-plane and along the 
c-axis. 

a. the current in the a6-plane 

The gap equation for the order parameter in the pres- 
ence of a uniform current can be written as 

A(k) = -T^^y(k,p)Tr[piaiG(ic.„,p)] (5) 

where A(k) = A/(k), V{k,p) = l//(k)/(p)*, and 
/(k) represents the momentum dependence of the or- 
der parameter. Here the single particle Green's function 
G'(ia'„,k) is given by 

G"^(zcj„,k) = iujn - Vk • qs + CkPa + A(k)picri, (6) 

where ^k = (fc^ + fc^)/2m — tcos{ckz) — ^ is the single 
particle dispersion. 



The uniform current reduces the amplitude of the order 
parameter and the current dependence of the amplitude 
at r = can be obtained from the following equation 
derived from Eq.4 (the lattice constant c = 1 for simplic- 
ity)- 



In 



A(0) 
A(.) 



^ / dk 



d(h I f P Re arccosh , , , , 

A(s)|/|_ 

(7) 

where s — vpqs, is the angle between the direction 
of the current and the quasiparticle velocity Vk. A(s) 
and A(0) represent the amplitude of the order parameter 
in the presence and absence of the current, respectively. 






0.4 
S/A(0) 



FIG. 1: The supercurrent as a function of s/A(0) with the 
current (a) in the afe-plane for both of the /-wave and /'-wave, 
(b) along the c-axis for the /-wave, (c) along the c-axis for the 
/'-wave. The unit of the supercurrent is {en/mvF)A{0) and 
{en/m){vFc/vp)A{0) for the current in the plane and along 
the c-axis, respectively. 



In the case of the /-wave order parameter, /(k) = leads to 
e^*'''' cos(cfcz), a straightforward computation at T — 



hi 



rA(o)i 
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[a(s)J 


TT 



arcsinj/ \ny 



arcsm y 



d(j) ln(sin( 



(8) 



for y = s/A(s) < 1. When y > 1, there is no solution for the gap equation. Now the contribution from the 
quasiparticles to the current can be computed from '33| 



hp = eT^^Tr[vkG'(iw„,k)]. 



(9) 



Taking into account this, the net supercurrent for the f-wave order parameter is obtained as 



en 
m 
en 
m 



dk. 



T^^V Jo 



1 



Try 



d(f) COS0 Re\/(ycos0)2 — (cosfc^)^ 

l-2y^ . 
arcsm y 



(10) 



where n — kp /2Tr is the density of electrons in the plane. The supercurrent as a function of s/A(0) can be obtained 
from EqlSjand EallOl the result is plotted in Fig. 1(a). The critical current occurs at s = 0.74A(0) and the value of 
the critical current is jc — 0.465(en/mtii?) A(0). 
b. the current along the c-axis 

When the current is along the c-axis, the gap equation for the /-wave order parameter at T = is now given by 



In 



rA(o)i 




[a(,s)J 


TT Jo 



dkz (cos/cz)^ Re arccosh 



s sin kz 
A cos kz 



— arcsinh y - 



where s = vpcQs with vpc — tc and Vk = (te) sinfe^q^ is used. The net supercurrent is also found as 



en / vpc 



m \ vp 



1 / dkz sin kz Re \/ ( y sin kzV ~ (cos kzV 

T^y Jo 



en I vpc 



1 - 



(11) 



(12) 



Notice that vpc and vp are the velocities along the c-direction and in the a5-plane, respectively. The supercurrent 
computed from Eg II II and Eq^jis plotted in Fig. 1(b); the critical current occurs at s = 0.61A(0) and the value of 
the critical current is jc — 0.277{en/ni){vpc/vp'^)A{0). 

/'-wave order parameter: When the current is in the a5-plane, the gap equation and the expression for the 
supercurrent turns out to be the same as those of the /-wave case. Thus the supercurrent in this case is given by Fig. 
1(a). On the other hand, when the current is along the c-axis, both of the gap equation and the supercurrent for the 
/'-wave order parameter are quite different from the /-wave counterpart. The gap equation now has the following 
form with s = vpcQs- 



In 



A(0) 



dk. 



ln(2y) - - I 1 



2 
1 



Re arccosh 
1 



s sin kz 
Acos(^) 



9[y- 



(13) 



where 9{x) — 1 for x > and is zero if x < 0. The supercurrent is given by 



en / vpc 

is = — 

m \vp 



1 / dkz sin kz Re\ (y sin kz)^ 

T^V J a V 



4 



The supercurrent as a function of s/A(0) is plotted in Fig. 1(c). Notice that the quasiparticle contribution to the 
current does not enter for s/A < 1/2, so the supercurrent is proportional to in this regime. The critical current is 
found at s = 0.56A(0) and the value of the critical current is jc — 0.531(en/m)(-(;i?c/w_F^)A(0). 



Summary and Conclusion: We investigated the be- 
havior of the supercurrent and obtained the critical cur- 
rent in the cases of the /-wave and /'-wave order param- 
eters defined in Eq.2 and Eq.3. In the superconductor 
with the full gap, the quasiparticle contribution to the 
supercurrent would enter only when the current exceeds 
some value. If the order parameter has nodes, however, 
the quasiparticle contribution may affect the supercur- 
rent even in the small current limit. When the current 
flows in the a5-plane, the supercurrent in the /-wave and 
/'-wave cases is indeed affected by quasiparticles even for 
small current and it has the same form in both cases. On 
the other hand, if the current is along the c-axis, the su- 
percurrent in the /'-wave case behaves similarly to the 
case with the full gap; in contrast, the quasiparticle con- 
tribution still enters at small current in the /-wave case. 

Now the discussion about different models for the su- 
perconducting order parameter in Sr2Ru04 is in order. 

1) When the current is in the a6-plane, the supercur- 
rent in the multigap model would be mainly determined 
by the dominant p-wave component on the 7-band and 
its behavior is similar to the case of s-wave superconduc- 
tor; the small contribution from the /'-wave component 
is subdominant so that the supercurrent is not much af- 
fected by the quasiparticles at small current. Thus, the 
critical current in the multigap model would be much 
bigger than that of the /-wave model where the super- 
current acquires a significant quasiparticle contribution. 

2) If the current is along the c-axis, even the supercur- 
rent from the /'-wave component (as well as the p-wave 
component) in the multigap model does not acquire the 
quasiparticle contribution at small current and as a result 
the supercurrent is even less affected by the quasiparti- 
cles. On the other hand, the supercurrent in the /-wave 
model is still very much affected so that the critical cur- 
rent is again much smaller. 

The qualitative difference in the behavior of the super- 
current in two models can be used to discriminate one 
of the leading candidates for the order parameter of the 
superconducting phase in Sr2Ru04. 
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